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ABSTRACT 

Angular momentum is transported outwards through an accretion disc by magneto- 
hydrodynamical (MHD) turbulence thus allowing material to accrete on to the central 
object. The magneto-rotational instability (MRI) requires a minimum ionisation frac- 
tion to drive turbulence in a disc. The inner parts of the disc around a young stellar 
object are sufficiently hot to be thermally ionised. Further out, cosmic rays ionise the 
surface layers and a dead zone forms at the mid-plane where the disc is too cool for the 
MRI to operate. The surface density in the turbulent active layer is often assumed to 
be constant with radius because the cosmic rays penetrate a constant layer. However, 
if a critical magnetic Reynolds number, i?eM,crit, is used to determine the extent of the 
dead zone, the surface density in the layer generally increases with radius. For small 
critical magnetic Reynolds number of order 1 , the constant layer approximation may 
be a reasonable fit. However, MHD simulations suggest the critical magnetic Reynolds 
number may be much larger, of order 10*. Analytical fits for the surface density in the 
magnetic active layer show that Sm oc Re^^^-^^R^^'^T^ , at temperature T and radius 
R, are a good fit for higher critical magnetic Reynolds number. For the metallicity 
variation between our galaxy, the LMC and the SMC, there should be no significant 
difference in the extent of the dead zone. Observations suggest an increase in the life- 
time of the disc with decreasing metallicity that cannot be explained by the dead zone 
structure (ignoring possible differences in dust abundances). 
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1 INTRODUCTION 



Low-mass stars are thought to form from the free-fall 
collapse of a protostellar molecular cloud core to a 
protostar with a disc on a timescale of a few 10^ yr 
l|Shu. Adams fc Lizand Il987l ). Angular momentum trans- 
port in accretion discs is driven by turbulence thus allow- 
ing material to accrete on to the young star. The magneto- 
rotational instability (MRI) can dri ve turbulence if the gas i s 
well coupled to the magnetic field (|Balbus fc Hawlevlll99ll ). 
However , with a low i onisation fraction the M RI is sup- 
pressed (|Gammiell 19961 : iGammie "fc MenoullTggi ). The inner 
parts of a disc around a young stellar object are hot enough 
to be thermally ionised. However, further out, the disc be- 
comes layered with an MRI turbulent (active) layer at each 
surface and a dead zone at the midplane that is shielded 
from the ionizing radiation of cosmic ray s and X-rays from 
the star (e.g. iSano et al..,200Q : Matsumura fc Pudritzll2003h . 
In this work we concentrate on circumstellar discs, but note 



that dead zone formation is also fav ourable in circumplane- 
tary discs (|Martin fc LubowllioTlal lbl'). 

A frequently used assumption in calculations 
of discs with dead zones is that the surface den- 
sity of the MRI a ctive surface layer is constant 



with radius (e.g. lArmitage. Livio fc Pringl^ 



Zhu. Hartmann fc Gammid 2009a : TerguenJ 



Matsumura. Pudritz fc Thommed \200^ . The 



2001 : 



20081 : 



rays enter the disc surface and are attenuated exponentially 
with a stopping depth of surf ace density around lOOgcm"^ 
(|Umebavashi fc Nakandll98lh . However, this does not allow 
for recombination effects that may play a significant role in 
determining the ionisation of the disc. 

A more realistic way to determine the dead zone is using 
a magnetic Reynolds number. MHD simulations show that 
magnetic turbulence cannot be sustained if the magnetic 
Reynold s number is lo wer than some critical value Rem < 
i?eM,crit (Flemins 2000"). However, the critical value is uncer- 
tain. Simulations of Fleming (2000) that compute the non- 
linear outcome of instability suggest that -RcM.crit ~ 10^ 
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without a net magnetic flux through the disc. With a mag - 
neti c flux it could be of the order of 100. .WardQ (|l999l ') 
and iBalbus fc TerauemI l|200ll ) considered the linear stabil- 
ity of the discs and find that the effects of Hall electromo- 
tive forces are important and the value may be as small 
as i?eM,crit = 1- However, conditions for linear instability 
and for maintaining turb ulence are not necessarily the same 
l|Balbus fc Hawlevll2000l ). We consider a range of values for 
the critical magnetic Reynolds number from 1 to 10'*. 

The surface la yers of the disc may be ionised by cos- 
mic rays or X-rays. iMatsumura fc Pudrit j (|2003t ) flnd that 
cosmic rays determine the extent of the dead zone under 
typical conditions unless the X-ray energy is very high. We 
consider only cosmic ray ionisation and neglect X-rays. Elec- 
trons in the disc are removed through dissociative recombi- 
nation with molecular ions and at a slower rate with radia- 
tive recombination with heavy metal ions. The ionisation 
fraction, and hence the extent of the dead zone, is sensitive 
to the number of metal atoms because they quickly pick up 
the charges of molecular ions but slowly recombine with the 
electrons. 

The typical life time of a protostellar disc in the so- 
lar neighbourhood, where the me tallicity is Z m 0.02 i s 
of the order of 3 - 5Myr (e.g. [Hernandez et al.1 |2009| ). 
For sta rs near SN 1987A in the Large Magella nic Cloud 
(LMC) iDe Marchi. Panaeia fc Romaniellol (j2010h find from 
Ha fluxes that the mean lifetime is of the order of 13 Myr 
and the metallicity there is Z ~ 0.008. Note however, that 
these stars are likely to be slightly mo r e mas sive. Sim- 
ilarly, IPe Marchi. Panaeia fc Romaniellol ()201ll ) observed 
the NGC 346 cluster in the SmaU Magellanic Cloud (SMC), 
where the metallicity is Z ~ 0.002, and find disc lifetimes 
of around 20 Myr. It appears from these results that the life 
time of the disc increases with decreasing metallicity. Dead 
zones affect the accretion rate through the disc and thus 
lifetime of disc accretion and so we consider how the dead 
zone varies with metallicity. 

In Section [5] we describe the layered disc model where 
the extent of the dead zone is determined with the critical 
magnetic Reynolds number. We consider thermal ionisation, 
cosmic ray ionisation and recombination. In Section [3] we 
choose a fiducial disc model and find the dead zone and re- 
sulting accretion rate through the disc for varying metallic- 
ity and critical magnetic Reynolds number. We investigate 
the extent to which assuming a constant surface density in 
the active layer can approximate this and how sensitive the 
dead zone is to the metallicity. In Section 3] we find analyti- 
cal approximations to the surface density in the active layer 
of the disc. 



2 LAYERED DISC MODEL 



tions of MHD turbulence find a > 0.01 (^Brandenburg et~al 



In an accretion disc l|Lvnden-Bell fc Pringld Il974l : iPringlel 
Il98ll ) around a central object of mass M, where the material 
is in Keplerian orbits with angular velocity Q, = ^ GM/R^ 
at radius R, the kinematic turbulent viscosity may be 
parametrised by the Q-prescription 

-4 w 



jShakura fc Sunvaev|[l973 ). where q is a dimensionless pa- 
rameter that is not well constrained. Numerical simula- 



1995: Stone et al 



20091 : iDavis et al. 



.19961 : iFromang et~al. 2007; Gu an et al . 
2009h . Observational evidence from dwarf 
nova and X-ray outburst s suggests a ~ 0.1 — 0.4 
dKing. Pringle fc Lividl2007l ). In this work we take a = 0.1 
but discuss implications of a lower value. The sound speed 
is given by 



(2) 



where T is the midplane temperature of the disc, TZ is the 
gas constant, /i = 2.3 is the mean particle weight. The disc 
is viscous only in the turbulent active surface layer of the 
disc, not in the dead zone. 

The dead zone is determined by where the MRI is un- 
sustained because of poor coupling of the magnetic field to 
the disc. If turbulence is damped on a scale of A, then the 
growth rate of the MRI is Va/A, where Va = B/{A-Kp)2 
is the Alvon speed , B is the magn e tic fie ld and p is the 
density. .Hawley. Gammie fc BalbusI (|l995l ) use numerical 
simulations to find the shear stress due to turbulence in- 
duced by the instab i lity to be Wrcp ~ pV^. In terms of the 
IShakura fc SunvaevI l| 19731 ) a-parameter this is = ctpc^ 
and so equating these we find Va ~ a^Cs. The growth rate 
of the Ohmic diffusion is rj/}^, where 77 is the diffusivity 
of the magnetic field. Complete damping occurs when all 
turbulence scales less than the scale height, H — Cs/0., are 
damped. This means that X < H, or equivalently we can 
define the magnetic Reynolds number 



Re^ 



a^CsH 



(3) 



Note that 



and the disc is dead if ReuiR, z) < Rbm 
several definitions of th e magnet i c Rey nolds number have 
been used. For example, iFleming (|2003) defined it as above 
but without the factor of a 2 . We consider a range of values 
for -ReM.crit- The Ohmic resistivity is 



234 T 2 2 
T] — cm ! 



(|Blaes fc Balbusll 19941 ). where the electron fraction is 



(4) 



(5) 



where rio is the number density of electrons and Tin is the 
total number density. With vertical hydrostatic equilibrium 
the number density is 



n-n ~ Tic exp 



(6) 



where z is the height above the mid-plane of the disc. The 
number density at the midplane of the disc, where z = 0, is 

(7) 



y27vpmuH 

where mn is the mass of hydrogen. In the following two 
sections we consider the electron fraction in the disc with 
thermal ionisation, cosmic ray ionisation and recombination. 



2.1 Thermal Ionisation 

The dominant ther mal ionisation ions in protos tellar discs 
are Na^ and K"*" (lUmebavashi fc Nakand Il98ll ). However 
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the ion is more important at the onset of dynamically 
interesting ionisation levels because of its smaller ionisation 
potential. The Saha equation for the electron fraction from 
thermal ionisation can be approximated by 



Xc = 6.47 X 10" 



2.4 X 10^ 



10 



[K/H]\ 2 



T 



103 K 



exp(-25188/T) 
1.15 X 10-" 



(8) 



iBalbus fc HawlevI |2000| ). where [K/H] = log^o{K/H) - 
logj^g(isr/ii")sc)iar is the potassium abundance relative to hy- 
drogen and log]^g(-ftr/-ff)soiar = —7. We consider the effect of 
varying [K/ H] in Sections 13.11 Thermal ionisation is most 
important in the inner parts of the disc where the tempera- 
tures are higher and it falls off exponentially. 



2.2 Recombination 

When thermal ionisation becomes negligible, the ionisation 
fraction is found with the balance of cosmic ray ionisation 
and recombination effects. Electrons may be captured by 
dissociative recombination with molecular ions (of density 
nj^+) and radiative recombination with heavy-metal ions 
(of density n„-j+). Charge will be transferred from molecular 
ions to metal atoms and so the electron fraction depends 
also upon the metal fraction 



Tin 



(9) 



where is the number density of metals. The standard 
notation for metallicity is by mass fraction, Z. This is related 
to the metal fraction with 



Z 



(10) 



where fj. is the average mass of a particle in units of the 
mass of hydrogen. For example, the molecule of highest mass 
taking part in the reactions may be CO, which has a mass 
of 28 units. In our galaxy Z = 0.02 so 7 x lO"**. In the 

LMC Z = 0.008 is equivalent to > 3 x lO""* and in the 
SMC Z = 0.002 is equivalent to Xn^ > 7 x 10"^ 

Assuming the rates are the same for all species, the rate 
equation for the electron density is 



(^dWM+ + /3r?lm+)"e 



and for the molecular ion density 

rf('T-M+) 



dt 



(11) 



(12) 



The ionisation rate, i^, is discussed in the next section. The 
recombination rate coefficients are 



/3r = 3 X 10~"r~3 cm^s"^ 



(13) 



for the radiative recombination of electrons with metal ions, 

(14) 



;Sd = 3 X IO'^T'2 cm^g-i 



for the dissociative recombination of electrons with molecu- 
lar ions and 



/3t = 3 X lO'^'cm^s"^ 



(15) 



for the charge transfer from molecular ions to metal atoms. 
Conservation of charge tells us that 



(16) 



Solving equations (|ll|l and (|12p in steady state with equa- 
tion (|16|) gives the cubic equation 



„3 I ^ 2 







llQppenheimer fc Dalgariiol Il974l : ?; iMillaJ 

iFarauhar fc Willacvl l that can be solved to find the 
electron density. There are three solutions to equation (|17|l 
but only one that is both real and positive. 

The solution to equation (I17|l can be simplified for the 
case where there is zero metallicity, Xm = 0, and find 



(18) 



and also when the metals dominate so that Xm ^ Xc and 



(19) 



(e.g. iMatsumura fc Pudritd [20031 ). We use these limits in 
Section [31 We note that this model does not take the effects 
of dust in to account that may increas e the size of the de ad 
zone (e.g. [Turner "fc SanQ,200a : ,Okuzumi fc Hirosdiioill ). 

2.3 Cosmic Ray Ionisation 

We consider the only external ionisation source to be cosmic 
rays with ionisation rate 

En,/2^ 



C = Co exp 



Xcr 



(20) 



(where we have n eglected the seco nd term in the equation by 
ISano et al.|[200ol ). Co = 10"^^ s~^ llSpitzer fc Tomaskdligesh 
and Xcr = 100 gcm"^ l|Umebavashi fc Nakandll98lD . The 
total surface density in the active layers is 



2/imH 



nndz = Eerfc ^-^||:^ , (21) 



where Za-it is the height of the disc above the mid-plane 
where the dead zone ends which we find in Section [231 



2.4 X-ray Ionisation 

You ng stellar objects may be active X-ray sources 



Kovama et al 



biects : 



However, IMatsumura fc Pudrit j 
(|2003f ) find that cosmic ray ionisation dominates X-ray ion- 
isation. For example, comparing their figures 2 and 3 the 
dead zone is much larger in both height and radius with 
only X-rays. X-rays may dominate only if the energy is very 
high, fcTx — 5— 10 keV, but this is much higher than for most 
observed sources. In this work we neglect X-ray ionisation 
and concentrate on the effects of cosmic ray ionisation. 



2.5 Dead Zone 

We have found the electron fraction with thermal ionisation 
and cosmic ray ionisation with recombination. In order to 
determine the dead zone height we take the electron fraction 
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to be the maximum of equation ((HJ and the solution to equa- 
tion (I17|) . Now we can find the magnetic Reynolds number in 
equation Q and for each radius, R, find the critical height 
in the disc, Zcrit by solving Rem = -ReM.crit- Where thermal 
ionisation is dominant, this is straightforward. However, it is 
slightly more complicated where cosmic ray ionisation dom- 
inates. For example, with the high metallicity limit for the 
electron fraction (equation I19|l we solve 



1 / ^crit 

2 KIT 



= ln 



(22) 



- ^erf (Ip^) . 

In the limit of low metallicity we replace /?r with /3d. Given 
a total surface density and temperature, this equation can 
be solved numerically to find Zcrit and hence the active layer 
surface density given in equation (pT|) and the dead zone 
surface density Ed = E — Em. We illustrate this with an 
example in the next section. 



3 DISC PROPERTIES 

We choose a fiducial model for the disc structure around and 
star of mass M — 1 Mq and consider properties of the dead 
zone. We take the total surface density to be 



(- 

VAU 



(23) 



where Eo = 10'^ g cm ^ and the midplane temperature is 



T : 



(24) 



where To — 100 K. We choose the disc such that without the 
dead zone, it has a constant accretion rate at all radii. This 
requires M oc i^E = const and so we take s + t = 1.5. Unless 
otherwise stated, we choose s — t — 0.75. This disc is in a 
steady state if it is supplied with a constant source of mass 
in the outer parts and the disc is fully MRI active. However, 
we now show that the disc may have a dead zone. 

3.1 Active Layer 

In Fig. [T] we show, at a fixed radius in the disc, R = 1 AU, the 
surface density in the active layer as a function of the total 
surface density. We vary the metallicity and critical mag- 
netic Reynolds number in the disc. The dotted lines show 
the high metallicity limit for the electron fraction given in 
equation (|19p . The solid lines show intermediate metallici- 
ties found by solving the cubic equation p7|l for the electron 
fraction. For fixed radius, the active layer reaches a constant 
surface density as the total surface density increases. The 
active layer is smallest with a low metallicity and a high 
critical magnetic Reynolds number. 

However, the finite metallicities considered in Fig. [1] are 
tiny. Even the highest finite metal fraction we consider here, 
Xm = lO"'^, corresponds to a very small metallicity of about 
Z — 3 X 10"" (see equation I10[) that is several orders of 
magnitude smaller than that in the SMC. The metallicity 
required for the dead zone to be significantly different from 
the high metallicity limit is very small, unless the critical 



o 



O 




3 3.5 

log,o(i;/g cm-=) 

Figure 1. For a fixed radius in the disc of i? = 1 AU, the surface 
density of the active layer in the disc is shown as a function of 
the total surface density. The upper two lines show -Rejvi.crit = !> 
the middle two lines show -Rem ^rit = 100 and the lower two 
i?eM crit = 10*. The dotted lines show the high metallicity limit 
for the electron fraction found in equation II19I I and the solid 
lines have Xm = 10~^^ (upper solid line), Xm = 10"^" (middle 
solid lino) and Xm = (lower solid line) found by solving the 

cubic H17|l for the electron fraction. 



magnetic Reynolds number is found to be significantly larger 
than 10*. Hence, for reasonable metallicities, those in the ob- 
servable universe, the dead zone structure does not depend 
on metallicity and the high metalliticy limit is appropriate. 
In the rest of this work we consider only the high metallicity 
limit. 

In Fig. [2] we show the height of the dead zone as a 
function of radius in the disc for varying critical magnetic 
Reynolds number in the high metallicity limit. As the crit- 
ical magnetic Reynolds number increases, the size of the 
dead zone increases. The re sults are similar to those of 
iMatsumura fc Pudritj liooS) who find the dead zone in a 
disc with a different surface density and temperature distri- 
bution. In Fig. |3] we also plot the inner most parts of the 
dead zone and see that the effect of varying the constant 
[K/H] (the Potassium abundance relative to hydrogen, see 
Section I2.1|l is not very significant. It only affects the in- 
nermost parts of the disc where thermal ionisation is the 
dominant source. Further out where cosmic ray ionisation 
becomes dominant the dead zone does not depend on [K/H]. 
In the rest of the work we take the solar abundance value, 
[K/H] = 0. 

In Fig.|4]we show the surface density of the active layer 
as a function of radius of the disc for varying critical mag- 
netic Reynolds number in the high metallicity limit. The 
inner parts of the disc are thermally ionised and the outer 
parts are fully ionised by cosmic rays. Hence in these re- 
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R/AU 

Figure 2. The height of the dead zone as a function of radius for 
critical magnetic Reynolds number of 1 (dotted line), 100 (short- 
dashed line) and lO** (long-dashed line) in the high metallicity 
limit. The solid line shows the scale height of the disc, H. Note 
that the axis scales are not the same. 

gions, the surface density of the active layer is equal to the 
total surface density, Em = E. 

In the region where a dead zone exists, the surface den- 
sity in the active layer initially decreases sharply. This is 
because the temperature in the disc decreases with radius 
(T cx i?""^''^) and the thermal ionisation decreases expo- 
nentially with temperature. However, once the ionisation 
through cosmic rays dominates the thermal ionisation, the 
active layer surface density increases with radius until the 
whole disc is turbulent. Independently, the decreasing tem- 
perature distribution has the effect of lowering the active 
layer surface density in this region. However, the dominant 
effect is the scale height of the disc, H, that increases with 
radius (H oc R^'"^). As described in Section[2l the turbulence 
is completely damped when it is damped on all scales < H . 
The damping decreases with radius because of the increas- 
ing scale height and thus the active layer surface density 
increases. 

The constant surface density in the active layer that is 
more generally assumed is also plotted in the dot-dashed 
line in Fig. |4l We take as an example a disc which is dead 
if r < r„it = 1500 K and E > E„it = 200 gcm'^. The 
active surface density determined by the critical magnetic 
Reynolds number is very different from a constant except for 
a the low critical magnetic Reynolds number. For i?eM,crit = 
1, the surface density of the active layer is close to constant 
but much larger than the generally assumed values that are 
less than around 200gcm~^. The minimum surface density 
in the active layer, that occurs close to the inner edge of 
the dead zone, can be very small especially for high critical 
magnetic Reynolds number. This is not well represented by 
the assumed constant value. 




R/AU 

Figure 3. The height of the inner parts of the dead zone as 
a function of radius for a critical magnetic Reynolds number of 
100 in each case with solar abundance, [K/H] = 0, (solid line, 
as used for the rest of the work), [K/H] = 2 (dotted line) and 
[K/H] = —2 (dashed line) in the high metallicity limit. The three 
dead zone heights are different where the electron fraction due to 
thermal ionisation (equation [Sj is larger than that due to the 
cosmic ray ionisation (equation 11911 . At larger radii thermal ion- 
isation becomes negligible and the dead zone is found with the 
cosmic ray ionisation that does not depend on the parameter a. 
Note that the scales of the axes are not the same. 

3.2 Accretion Rate 

The accretion rate through the active layer of the disc is 
approximately given by 

M = STTi/Em. (25) 

In Fig. Owe show the accretion rate as a function of radius in 
the disc for our fiducial model for varying critical magnetic 
Reynolds number in the high metallicity limit. The thin ac- 
tive layer severely limits the accretion rate over the dead 
zone by up to several orders of magnitude compared with 
that through the outer turbulent parts of the disc. These 
discs cannot be in steady state. 

The mass in the dead zone must increase in time be- 
cause it is continually supplied by the active layer. With 
sufficient accretion on to the disc, the outer parts of the 
dead zone can become gravitationally unstable. The extra 
heating from the viscosity due to self gravity can trigger the 
MRI and with the sudden increase in the turbulence a signif- 
icant fraction of the disc is accreted in an outburst. After the 
outburs t the disc cools, the dead zone reforms and the cycl e 
repeats (|Armitage, Livio fc Pringle|[200ll : IZhu et al.|[2009bl) . 
This outburst mechanism, known as the gravo-magneto in- 
stability, should be similar no matter how the extent of the 
dead zone is determined. The outbursts can be explained 
by plotting the accretion rate through the disc against the 
surface density at a fixed radius. There are two possible 



© 2011 RAS, MNRAS 000, [TH9] 



6 R. G. Martin et al. 




log,„(R/AU) log,„(R/AU) 



Figure 4. The total surface density (solid line) and the surface 
density in the active layer as a function of radius in the disc for 
critical magnetic Reynolds number of 1 (dotted line), 100 (short- 
dashed line) and 10* (long-dashed line) in the high metallicity 
limit. The inner parts of the disc are thermally ionised so Zjm — Zj 
and similarly the outer parts of the disc are fully ionised by cosmic 
rays. The dot-dashed line shows the surface density in the active 
layer for the constant layer assumption where the disc is dead if 
T < Tcrit = 1500 K and S > Ecrit = 200 gcm-^. 



Figure 5. The accretion rate through the active layer of the disc 
in the high metallicity limit. The solid line shows a disc that is 
fully turbulent without a dead zone. The other three lines have a 
dead zone with a critical magnetic Reynolds number of 1 (dotted 
line), 100 (short-dashed line) and 10* (long-dashed line). These 
discs cannot be in steady state. The dot-dashed line shows the 
accretion rate through a disc with the constant surface density in 
the active layer (as in Fig.[4]|. 



steady disc solutions, one that is fully turbulent and a sec- 
ond that is self gravitating. Outbursts occur when there is 
no steady solution for a given accretion rate at a specific 
radius . The disc m oves between the two solutions in a limit 
cycle l|Martin fc L ubow 201l3). 

The typical bolometric luminosity of protostars is 
smaller than would be expected fro m the infall rate an d 
duration of the protostellar phase (jKenvon et al.l Il990l ). 
Observations of T Tauri stars suggest that at an age of 
around 1 Myr the ac cretion rate th r ough t he disc is around 
10"^ Me yr~^ (e.g. IValenti et all 1 19931 : iHartmann et all 
Il998l ). One solution to this luminosity problem is thought 
to be time-dependent accretion on to the star that may be 
caused by dead zone formation in the disc. The disc spends 
the majority of its time with a dead zone and a restricted ac- 
cretion rate on to the star. A low critical magnetic Reynolds 
number of around one could help to explain the observed ac- 
cretion rates (see Fig. [SJ. In a time dependent disc, the ac- 
cretion on to the central star will be of the order of the 
minimum accretion rate through the disc after a viscous 
timescale. 



V 



(26) 



at that radius. For example, at a radius of 0.1 AU, the vis- 
cous timescale is short at around 450 yr. It is difficult to 
reconcile a model with a high critical magnetic Reynolds 



number to the observations of T Tauri star accretion rates 
unless there is a source of external heating on the disc. 



4 ANALYTIC APPROXIMATIONS 

Fig. m suggests that the surface density in the active layer 
above the dead zone can be approximated by a power law in 
radius for high critical magnetic Reynolds numbers. When 
the surface density is sufficiently large, the active layer sur- 
face density is insensitive to the total surface density (as 
shown in Fig. [l|. This, however, does not translate into the 
active layer surface density being constant in radius. As de- 
scribed in Section FS.ll the dominant effect is the scale height 
of the disc, H, that increases with radius leading to an in- 
crease in the active layer surface density. 

An analytic approximation to the active layer surface 
density is found by neglecting the third term in equation (|22l) 
and then approximating Em oc exp[— (zcrit/-ff)^/2]. For cos- 
mic ray ionisation 



1.36 



0.1/ Viooy V 10* 

lAuJ 



M 



Co 



g cm 



(27) 



In regions where thermal ionisation is the dominant source. 
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log,„(R/AU) 



Figure 6. The surface density in tlie active layer. The solid lines 
show numerical solutions. The dashed lines show analytic approx- 
imations. The upper lines have -ReM.crit = 100 and the lower lines 
have -Rejvi.crit = 10^. 



the layer can be approximated with 



- 0-2 (^) 



T ^ / -ReM.crit 
104 



M 



1500 
exp 



R 



0.03 AU 



9/2 



-25188 - 



r 1500 



g cm 



(28) 

Now we can fit the whole active layer surface density ana- 
lytically. In Fig. |6] we show the analytic approximation to 
the surface density in the active layer is a good fit for high 
critical magnetic Reynolds number, 7?eM,crit > 100. Note 
that even though we have only shown one temperature and 
surface density distribution in this work, these fits are valid 
for any distribution. 

The accretion rate through the layer can be approxi- 
mated by 
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cosmic ray ionisation is the dominant 
source. The minimum accretion rate through the disc occurs 
where cosmic ray ionisation takes over from thermal ionisa- 
tion. By finding where the surface densities in equations \27\ 
and (|28[l are equal, we find this occurs at a temperature of 
T ~ 10^ K. The radius corresponding to this temperature 
has the smallest accretion rate (see Fig. [5]). As explained in 
Section 13.21 the accretion on to the central object will be 
limited by this minimum over a viscous timescale at that 
radius. 

If the viscosity a parameter is smaller than the value we 



have taken in this work, of 0.1, then both the active layer 
surface density and accretion rate would be smaller than 
predicted here (see equations 1271 and I28p . The disc would 
become even more unstable to the gravo-magneto instabil- 
ity. The results shown here represent an upper limit to the 
surface density in the active layer. 



5 DISCUSSION 

There are three non-ideal MHD effects in the generalised 
Ohm's law, Ohmic resistivity, Hall effect and ambipolar 
diffusion, that dominate for different regimes. The Ohmic 
term, that we have used here, dominates at high density and 
very low ionisation. However, the ambipolar diffusion domi- 
nates for low density and high ionisation and the H all effect 
dominates in a region between these two extremes (|Wardlel 
Il997h . Recent work suggests that the MRI in the sur- 
face layer s of discs may be d ictated by ambipolar diffusion 
( Perez-Becker fc Chian3l2011a.b; Bai fc Stone 2011.) or the 
Hall effect (|Wardle fc SalmeronI l201ll ) rather than Ohmic 
resistivity as we have assumed here. IWardle &: SalmeronI 
(2011) applied a model with all three effects to the minimum 
mass solar nebula and found that at a radius of 1 AU there 
is around an order of magnitude increase or decrease in the 
active layer surface density, depending on whether the field 
is parallel or antiparallel to the rotation axis, respectively, 
compared with a disc with only Ohmic resistivity. These 
extra effects could change the surface density in the active 
layer and the accretion rate through the disc predicted here 
by around an order of magnitude. However, Ohmic resistiv- 
ity has the advantage of being independent of the magnetic 
field (Fleming 2000,) and the values it produces for the sur- 
face density of the active layer represent an average for the 
values obtained over a ran ge of vertical magnetic fields (e.g. 
IWardle fc Salmeronll201ll'). 



Fromang. Terguem fc Balbuj (|2002l ) use a similar ap- 



proach to the work presented here but find that the ac- 
tive layer in the disc decreases with radius such that the 
dead zone surface density is approximately constant. The 
temperature distribution they use, shown in their Figure 1, 
has a very steep radial dependence in the region where the 
dead zone forms. With the analytical fits given here, it is 
clear that the constant active layer surface density can be 
recovered if the temperature of the disc has a radius depen- 
dence T oc J?~®'^4. If the power is even smaller, < —9/4, 
then the active layer will dec rease with radius as found by 
iFromang. Terguem fc Balbusi (,2002 ) . The analytical fits pre- 
sented here in Section |4] would provide an approximation to 
the disc in this case too because they can describe any sur- 
face density and temperature distribution. 

If cosmic rays are present, the depth of the active layer 
is insensitive to the presence of X-rays because they do not 
penetrate so far. However, X-rays could still increase the 
ionisation fractio n, and h ence the associated effective a, in 
the active layer. Ilgea fc Glassgold (1999) showed that the 
X-ray ionisation rate near the surface is significantly larger 
than that of cosmic rays. The accretion rate through the 
layer may be higher than predicted in equation (|29|l if this 
is taken into account and perhaps could help to explain the 
T Tauri accretion rates for higher critical magnetic Reynolds 
number. 
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Observed FU Orionis systems may be one manifesta- 
tion of the time-dependent accretion phenomenon. These 
luminou s young stellar objects are found in star forming 
regions l|Hartmann fc KenvonI IT996I ) . Their optical bright- 
ness increases by five magnitudes or more on a ti mescale 
of ye ars and decays on a timescale of 50-100 years (|Herbiel 
Il977h . The timescale and brightness of accretion outbursts 
from the gravo-magneto instabi lity are similar to the ob- 
served FU Orionis ou tbursts (|Armitage. Livio k, Pringl3 
I2OOII : IZhu et al.l l2009bh . Time-dependent numerical simu- 
lations are needed to investigate these effects more fully in 
a disc with a dead zone determined by the critical mag- 
netic Reynolds number. By comparing the resulting out- 
bursts with FU Orionis observations the value of the critical 
magnetic Reynolds number may be further constrained. 



6 CONCLUSIONS 

A typical assumption used in time-dependent simulations of 
accretion discs with dead zones is that the surface density in 
the MRI active layer above the dead zone is constant with 
radius. However, when the dead zone is identified by the 
value of the critical magnetic Reynolds number including the 
effects of recombination, the active layer surface density is 
generally found to increase with radius. The constant layer 
is best reproduced with a low critical magnetic Reynolds 
number. 

However, MHD simulations suggest the critical mag- 
netic Reynolds number may be much higher. For higher 
critical magnetic Reynolds numbers, i?eM,crit ^ 100 we have 
found an analytical fit to the active surface density (equa- 
tions [27] and I28p that will be a useful approximation in fu- 
ture time-dependent calculations. The dead zone structure is 
very sensitive to the value of the critical magnetic Reynolds 
number, as seen in Fig. U but the value is still uncertain and 
needs to be clarified in further work. However, this is com- 
plicated by the fact that the a parameter in the viscosity is 
also uncertain. 

The metallicity variation between our galaxy, the LMC 
and the SMC is not significant enough to affect the expected 
size of the dead zone in a disc (ignoring possible differences 
in dust abundances). When comparing accretion discs in our 
galaxy with those in the LMC and SMC for example, the 
high metallicity limit will be appropriate. In order to explain 
the observed difference in disc lifetimes w i th me tallicity (e.g. 
iDe Marchi. Panagia fc RomanieUcl |2010| . I2OIII ') some other 
effect must be taken into account. 
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